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ABSTRACT 



Spitzer images of Tycho's supernova remnant in the mid-infrared reveal limb- 
brightened emission from the entire periphery of the shell and faint filamentary struc- 
tures in the interior. As with other young remnants, this emission is produced by dust 
grains, warmed to ~ 100 K in the post-shock environment by collisions with ener- 
getic electrons and ions. The ratio of the 70 to 24 /xm fluxes is a diagnostic of the 
dust temperature, which in turn is a sensitive function of the plasma density. We find 
significant variations in the 70/24 flux ratio around the periphery of Tycho's forward 
shock, implying order-of-magnitude variations in density. While some of these are 
likely localized interactions with dense clumps of the interstellar medium, we find an 
overall gradient in the ambient density surrounding Tycho, with densities 3-10 times 
higher in the NE than in the SW. This large density gradient is qualitatively consis- 
tent with the variations in the proper motion of the shock observed in radio and X-ray 
studies. Overall, the mean ISM density around Tycho is quite low (~ 0.1-0.2 cm"^), 
consistent with the lack of thermal X-ray emission observed at the forward shock. We 
perform two-dimensional hydrodynamic simulations of a Type la SN expanding into a 
density gradient in the ISM, and find that the overall round shape of the remnant is still 
easily acheivable, even for explosions into significant gradients. However, this leads 
to an offset of the center of the explosion from the geometric center of the remnant of 
up to 20%, although lower values of 10% are preferred. The best match with hydro- 
dynamical simulations is achieved if Tycho is located at a large (3-4 kpc) distance in 
a medium with a mean preshock density of ~ 0.2 cm"^. Such preshock densities are 
obtained for highly (> 50%) porous ISM grains. 
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1. Introduction 



Tycho's supernova remnant (SNR; hereafter Tycho), also known as G120.1+1.4, 3C10, 
and Cassiopeia B, is a member of a small subclass of Galactic SNRs k nown as the "historica l 



supem ovae" (SNe). Tycho is the remnant of the SN observed in 15 72 dStephenson & Green 



20021), and was first characterized as a "Type I" SN 



ex plosion were first dis covered by 



Rest et al, 



by 



Baade 



1945|) . Light echoes from Tycho's 



by 



Krause et al. 



teOOSh . and subsequent spectroscopy of the echoes 



(120081) determined that Tycho resulted from a standard Type la SN. Various 



distances have been proposed in the lit erature for Tycho; an extensive comparison of various 



distance determinations is present ed in 



distance of 2.3 kpc, suggested by 



Havato et al. 



Chevalier et al.l ( 



mm 



distances in excess of 4 kpc have also been reported (jSchwarz et al 



Many authors have ad opted the 



1980) and 



Albinson et al.l (|l986|) : however, 



19951). We adopt a distance of 



2.3 kpc for this work, but also examine the effects of larger distances. 

The remnant has been widely studied across the electro magnetic spectrum . At radio 



wavelengths, Tycho exhibits a classic shell-type morphology (IDickel et al. 



19821) . while at optical 



wavelengths, th e shock is Balmer-dominated, with emission only seen in Ha from the eastern a nd 



northern limbs (jKamper & van den Bergh 



19781: 



Ghavamian et al. 



2000). 



Ishihara et al. 



( 20101) 



observed the remnant with AKARI between 9-140 /im, attributing the emission to dust in the ISM 
swept up by the forward shock. They report the possibilit y of as many as a few tenths of a solar 



mass in dust. In the far-infrared (IR), 



Gomez et al. 



(|2012h reported a clear detection of the entire 



shell at 70 and 100 fim with Herschel, finding the integrated flux from the remnant to be consistent 
with a small amount (< IO'^Mq) of warm dust at ~ 90 K. No clear emission associated with 
Tycho was observed at 160 /xm or beyond. In X-rays, a more complicated pi cture emerges, with 



the forward shock marked by thin rims of nonthermal synchrotron emission (Hwang et al. 



2002) 



he X-ray emission from the interior of the remnant is dominated by strong lines of Fe, Si, and S 



dWarren et al 



2005 



Havato et al 



2010|) from the ejecta. However, thermal X-ray emission from 
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the for ward-shocked interstellar medium (ISM) material is conspicuously lacking. iHwang et al 



(I2OO2I) searched within the Chandra data for this thermal component, only finding it in a few 



small sections of the remnant, and placing upper limits elsewhere. Recently , the remnant was 



detected in 7-rays, at both 



with Fermi (Giordano et al. 



eV e nergies with VERITAS (lAcciari et al 



2012 ) 



201l|) and GeV energies 



Tycho has generally been considered a prototypical Type la SN, expanding into uniform 
surroundings. However the ex act density of these surroundings has been a matter of debate. 



Dwarkadas & Chevalierl (|l998|) examined Tycho in the context of one-dimensional (ID) 
explosion models of a white dwarf with an exponential ejecta density profile. They found that 
to approximately match the observed size and ex pansion rate, a n ambi ent density of 0.6-1.1 



Badenes et al. 



(1200 6l) . who compare the 



cm"^ is required. A similar result was found by 
spatially-integrated X-ray spectrum of the western half of the remnant with synthetic spectra 
produced by ID delayed detonation hydrodynamic models, finding a pre-shock density of ~ 1 



cm 



On the other hand, an argument for a lower density in the ambient medium is made by 



Cassam-Chena'i et al. 



(120071) . who examine the lack of thermal X-ray emission from the blast 
wave, concluding that the azimuthally-averaged density surrounding Tycho must be < 0.2-0.3 
cm"^ (although they note that systematic erro rs could push this lim it as high as 0.6 cm""*). This 



density estimate is consistent with that from 



Kirshner et al. 



(|1987|) . who detected broad and 



narrow H a emission arising from "knot g" along the eastern limb, finding a pre-shock density of 



0.3 cm-^ 



Dwarkadas 



(l2000h extended the modeli ng of Type la explosions with exponential ejecta 



profiles from ID to 2D, and 



Katsuda et al. 



(120 10|) used these results, in combination with a study 



of the X-ray p roper motions ov er a baseline of seven years, to place a limit on the density of < 0.2 



cm ^. Finally, 



Lee et al. 



(|2004) suggest, based on radio observations of a molecular cloud in the 



vicinity of the northeastern portion of Tycho, that the remnant may be expanding into a density 
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gradient (though they do not determine numerical values for the density). 

A departure fr om spherical symmetr y is also suggested by measurements of the proper 



motion of the shell. 



Revnoso et al. 



(119971) found azimuthal variations of a factor of three in the 



proper motions at 1 .4 GHz me asured around the pe ripherv of the shell, with those in the SW being 



higher l.han in the NE, though 



Moffett et al. 



(|2004|) reported some modification s to these resu 



HughesI (|2000|) confirmed these azimuthal variations from ROSAT images, and 



Katsuda et al 



ts. 



2010|) examined high-resolution Chandra images and det ermined that the X-ray p roper motions 



of the forward shock vary by about a factor of two. Also, .Ghavamian et al. 



(2000) stated that 



the eastern side of Tycho must be interacting with a warm ISM cloud because optical emission 
produced by photoionizing radiation from the forward shock could be seen tracing the outer edge 
of the cloud. Clearly, the ambient density surrounding Tycho is quite uncertain, and in this paper, 
we examine the ISM structure as implied by broad-band flux ratios in the IR. 

The Spitzer Space Telescope has provided a new window on the Universe in the mid- and 



far-IR. 5p;Yzgr has returned remarkab 



2007 



Hines et al 



20041 : 



Temim et al. 



e images of SNRs, particular ly at 24 and 70 /zm ([Blair et al 



20061 : 



Williams et al. 



201 Ibh . Young SNRs, like Tycho, 



are typically in the non-radiative portion of their evolution, meaning that the temperatures in the 
post-shock gas are high enough that radiative cooling via coUisionally-excited lines of metals 



(Williams et al. 


2012; 


Temim et al. 


2012 


), arising from warm dust, heatec 


in the post-shock gas 


via collisions with energetic ions and electrons ( 


Dwek 


1987; 


Dwek et al. 


1996; 


Williams et al. 



201 la() . In Type la SNRs, this dust is ambient ISM dust; no newly-formed ejecta dust has ever 
been found in these remnants. As we will show below, the temperature of the dust, and thus its 
relative emission at 70 and 24 /xm, is most sensitive to the density of the gas in the post-shock 
environment. Our goal for this paper is to use the Spitzer images of Tycho to directly measure the 
post-shock density at various azimuthal locations around the shell. Nonthermal radio emission 
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from young SNRs does not provide a direct measurement of the density, and thermal X-rays 
from the shocked ambient medium are virtually absent in Tycho (additionally, since thermal 
X-ray emission scales with the square of the density, estimates obtained from X-rays contain a 
dependence on the unknown filling fraction of the gas). Optical emission traces only the densest 
B aimer-dominated shocks in the NE and N. Thus, the IR observations, and subsequent modeling 
shown here, provide the best constraints yet on density variations in the ISM surrounding Tycho. 
Throughout this work, we quote the azimuthal angle as the position angle, east of north. 

2. Data 

Tycho was observed in December of 2004 with all three instruments on Spitzer (PI J. Rho, 
Program ID 3483). Here, we report only on the 24 and 70 /im images from the Multiband Imaging 
Photometer for Spitzer (MIPS). We use the Post-Basic Calibrated Data (PBCD) from the Spitzer 
pipeline, version 18.12, for our analysis. We show the 24 and 70 /im images in Figured] The 
entire limb-brightened shell is seen in both images, and additional filamentary structures that run 
roughly in a N-S direction are visible at 24 /xm. Significant variations in brightness are seen in 
the remnant as well, with regions in the northwest and east having surface brightnesses several 
times higher than elsewhere in the remnant. We discuss the flux variations further in the following 
sections. For further comparison, we show a high-resolution optical Ha image (courtesy of P.F. 
Winkler), which shows emission only from the eastern and northern limbs. 



-7- 



3. Results 
3.1. Measurements 

We divided the outer shell of the remnant into equally-sized segments for our analysis, 
where each section is 30" in the radial direction and 60" in the tangential direction. The radial 
dimension is chosen to be slightly larger than the point-spread function of Spitzer at 70 //m, which 
is approximately 20", to ensure that all flux from the rim is captured within our regions. We 
were able to fit 19 non-overlapping regions around the periphery of the shell, ensuring that the 
outer boundary of each region extends slightly beyond the extent of the IR emission in the radial 
direction. We show all 19 regions, plotted on top of a 3-color IR mosaic, in Figure |2l This mosaic 
contains, for visualization purposes only, the 12 /xm image from the Wide-Field Infrared Survey 
Explorer (WISE). 

To measure the IR flux ratio, we had to first decide which IR images to use. We had six 
choices: 12 and 22 /im from WISE, 24 and 70 /im from Spitzer, and 70 and 100 pm from Herschel. 
We chose the Spitzer data. Spitzer images at 24 and 70 fim show clear emission from the entire 
shell of the remnant. The 70 /im image suffers from regularly spaced "striping" artifacts related 
to the scan direction of the telescope during mosaic observations. This is a well-known issue wi th 



2010 ) 



the Spitzer 70 pm MIPS detector, and one that we have encountered before (|Sankrit et al. 
In that paper, we found that despite their unsightly appearance, the stripes contribute an overall 
uncertainty of only 5% to the fluxes measured at 70 fim. Here, a similar analysis of the variations 
found in "on-stripe" and "off-stripe" flux measurements from relatively uniform regions of the 
remnant produced a similar result. In all the regions we tested, most of the variations were of 
order 5%; the largest effect we found was 8%. Thus, we did not attempt to correct for the striping 
pattem; rather, we simply added an additional 8% uncertainty term to all measured fluxes at 70 
/im. The MIPS Instrument Handbook lists calibration uncertainties on fluxes of extended sources 
at 4% and 7% for the 24 and 70 fim detectors, respectively. With the additional 8% uncertainty to 



-8- 



the 70 fim fluxes, we conservatively assume aU flux ratios to have a 16% uncertainty. 

Emission from dust grains at 12 fim arises from the very smallest grains that undergo 
temperature fluctuations and which are particularly prone to destruction through sputtering by 
thermal ions. The WISE 22 /im image is essentially identical to Spitzer data at 24 //m, but the 
spatial resolution and sensitivity of WISE are both lower than Spitzer. For this reason, we do not 
use either of the WISE images in our analysis. 

The Herschel data present two issues. First, while the shell of the remnant is clea rly detected 



Gomez et al. 



m\m . Thus, our 



at 70 /im, the detection at 100 /im is rather weak (see Figure 8 of |j 
decision is between the 70 /im data from Herschel and Spitzer. Ideally, we would use the Herschel 
data, since the spatial resolution of Herschel is several times better than Spitzer, making the 
Herschel 70 //m image of comparable spatial resolution to the Spitzer 24 u m image. However, th e 



Anianoetal 



( 2012 ) 



second issue with Herschel data is calibration uncertainty. Recent work by 
has revealed significant issues with the extended source photometry calibration of the PACS 
instrument on Herschel, which contains the 70 /tm camera. Because we are only doing two-point 
photometry to constrain the temperature of dust grains, it is crucial that the two fluxes that we use 
are calibrated as closely to each other as possible. 

With our 19 spatial regions defined, measuring the flux is straightforward. We first convolve 
th e 24 ixm image to th e resolution of the 70 /xm image using the convolution kernels provided 



m 



Gordon et al 



(|2008|) . We define a background that consists of four separate regions, each a 
few arcminutes outside of the remnant, to the NE, NW, SW, and SE of the shell. Although the 
background is fairly uniform in the immediate surroundings of Tycho, we average these four 
regions to create a single off-source background that we subtract from each flux measurement 
(scaled to the size of our extraction regions). 

We report our flux measurements at 24 and 70 //m in Table [U The 24 jim flux shows 
significant variations as a function of azimuthal angle around the periphery of the shell, being as 



-9- 



much as 15 times higher in the NW than in the SW. The 70 flux varies by only a factor of two 
around the shell. The different behavior of the two fluxes is quite significant: the remnant is not 
only brighter in some places than in others, but the ratio of the 70 fimflux to the 24 fimflux, and 
thus the temperature of the dust, varies by nearly an order of magnitude from one place to another 
in the remnant. It is the ratio of the 70 to 24 (xm flux that we fit with our models, described below, 
to determine the gas density behind the shock. 



3.2. Modeling 



As previously stated, the ratio of IR fluxes resulting from emission from warm dust grains is 



a diagnostic of the conditions of the X-ray emittin g plasma. We 



l ave dev eloped spectral models 



Williams et al. 



(l2011al) for a more complete 



for dust emission in SNRs; we refer the reader to 
description. Briefly, the spectrum emitted by a dust grain immersed in a hot plasma depends on 
the temperature and density of both the electrons and ions. The grain is heated by collisions with 
particles, with proto n and alpha particle collisions also slowly eroding the grain via sputtering 



(Nozawa et al. 



2006|) . These grains exist in the ISM encountered by the forward shock wave of the 



SNR, and are not newly-formed grains from the SN ejecta. Grain properties are important as well: 
large grains are heated to lower temperatures than small grains. The smallest grains (below a few 
nm in size) emit radiat ion quickly en ough that they cool back down to their ambient temperatures 



before being hit again dDrame 



2003|) . and these stochastic temperature fluctuations must be taken 
into account. Larger grains reach an equilibrium temperature determined by the balance between 
the coUisional grain heating rate and the grain radiative cooling rate. Smaller grains are more 
quickly destroyed via sputtering. The optical properties of the grain are also impor tant, e.g., the 
IR sp ectrum from a carbonaceous grain looks different from that of a silicate grain (jPraine & Lee 
19841) . 



Our models take all of the above effects into account. We use the dust grain size distributions 



- 10- 



of 



Weingartner & Draine 



(|2001|) . which contain a mix of carbonaceous and silicate grains ranging 
in size from 1 nm to 1 fim, in proportions appropriate for the ISM of the Milky Way. We use 
their favored model for the ISM in the Milky Way, with Ry (the ratio of visual extinction to 
reddening) of 3.1 and be (a parameterization of the carbon abundance in small grains) of 6.0, 
though we note that th e choice of grain-size distribution does not significantly affect the derived 



fits to the IR emission (Williams et al. 



20061). Once the grain properties are defined, what remains 



is to determine the plasma conditions. Because grains are heated by the forward- shocked gas 
of roughly cosmic abundances and we assume that the contribution of heavier elements to the 
gas density is small, we can assume that the post-shock electron density is 20% higher than the 
proton density, i.e. = l.lup, while the alpha particle density is one-tenth the proton density, i.e. 

Ua = O.lUp. 



3.2.1. Ion Temperature 

Determining the plasma temperatures for the various regions in Tycho is not necessarily 
straightforward, but we use the following approach. First, for the proton temperature, we assume 
the standard shock jump conditions derived from the Rankine-Hugoniot conditions: 



kTp = —mpV,, 



(1) 



where nip is the mass of a proton and is the sho ck velocity. To deter mine the sh ock velocity, w e 



use th e proper motion measurements from radio ([Reynoso et al 



1997h and X-ray (IKatsuda et al 



20101) studies, scaling the expansion rate to a distance of 2.3 kpc. Both of these studies report the 
proper motion of the forward shock as a function of azimuthal angle for the entire periphery of the 
remnant, and the agreement between the two is generally within 20% (although, in select locations, 
the discrepancy between the two measurements is as high as 40%). Because of this, we simply 
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use an average of the radio and X-ray values to determine the shock velocities. For the regions in 
which the X-ray and radio proper motions agree to within 10%, we assign an uncertainty of 10% 
to the shock velocity, approximately equal to the errors reported in both papers. In areas where the 
discrepancy is greater than this, we use the absolute values of the difference between the average 
value and that from each study for the uncertainty. The largest uncertainty in the shock velocity 
obtained in this way is 20%. We list all shock velocities, with uncertainties, in Table [T] along with 
the proton temperatures derived from these shock velocities. We assume that the downstream 



energy loss of the protons to Coulomb collisions has been minima 
for temperature equilibration found in the Ha-emitting shocks by 



, consistent with the low values 



Ghavamian et al 



(l2000f) . We 



also assume that alpha particle temperatures are four times higher than proton temperatures, i.e.. 



that ion-ion equilibration is negligible. This is consistent with results fro m SN 1006, w 



lere shocks 



of similar speed have been found to have minimal ion-ion equilibration (|Laming et al 



1996) 



There are several caveats to the proton temperature determination. First, as previously 
mentioned, temperature equilibration will bring the proton and electron temperatures closer 
together. However, after a fast initial rise in the electron temperature from ~ 10,000 K up to a 
few 100,000 K, the rem aining equilibration proceeds quite slowly, and is negligible in Tycho 



(IGhavamian et al. 



120011) . An additional effect is introduced by the uncertainty in the distance 
to Tycho. Lastly, the Rankine-Hugoniot conditions themselves assume that no energy is lost 
in the shock to escaping particles, an assumption which breaks down in the case of efficient 
co smic-ray acce l eration . Ty cho may well be the site of such acceleration, as has been proposed 



by 



Warren et al. 



(2005) and 



EriksenetaL 



(|2011|) . Efficient particle acceleration will lower the 
post-shock proton temperature from the values we report. A factor of tw o change in kTp s olely 



due to cosmic-ray acceleration would imply a compression ratio of ~ 10 (|Vink et al 



20101) . 



However, the effects of these uncertainties in the proton temperature on the modeling of 
warm dust emission are not large. Even a factor of two variation in kTp only results in a 25% 
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change in the inferred densities. We can quantify the dependence of the densities we report in 
Section |4] on the uncertainty in the distance. We assume D = 2.3 kpc, but if the distance were, for 
instance, 3 kpc, the proton temperatures implied by the shock velocities would be 70% higher. 
Higher temperatures would lead to lower densities, but the densities we report would only be 
lower by a factor of 1.2. Finally, it is important to note that these uncertainties in the proton 
temperature will affect only the absolute determinations of the density, and not the relative density 
from place to place in the remnant. 



3.2.2. Electron Temperature 



Hwang et al.l (120021) report electron temperatures of around 2 keV in Tycho based on a 49 ks 
Chandra observation in September of 2000. Atomic databases used in X-ray modeling codes like 
XSpec have been updated significantly since 2002. We reanalyze the X-ray data for this paper, 
using NEI version 2 .0, augmented by cus tom atomic line codes which include missing inner-shell 



electron transitions ( Badenes et al.l 



20061). Since Tycho was reobserved for s ignificantly longer 



in 2009, we use the new data, but repeat the analysis of 



Hwang et al. 



(2002). We also include 



XMM-Newton observations (PI A. Parmar) and fit both spectra independently. 

We derive a sl ightly lower tempera ture of 1.35 keV. To explore the discrepancy between our 



fits and those from 



Hwang et al. 



(|2002|) . we applied the spectral models used to fit the 2009 data 
to the original 2000 data, finding the same lower temperature. We obtain statistically identical fits 
from the XMM -Newton data. We believe the most likely explanation for the discrepancy between 



Hwang et al. 



(I2002|) and our work is simply a change in the calibration of the telescope and/or 
updated atomic data over the last decade. We therefore set the electron temperature at 1.35 keV, 
and make the assumption that this is the temperature in every region around the periphery of the 



shock. 
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Obviously, this assumption is unlikely to be correct everywhere along the rim. But how much 
does the value of the electron temperature matter for our dust modeling? As it tums out, it matters 
even less than the proton temperature does. A factor of two difference in the electron temperature 
has only a 10% effect on the ratio of the IR fl uxes measured in the Spitzer bands. The difference 



between 1.35 and 2 keV, the value reported in 



(l2002h . has only a 5% effect on the 



Hwang et al. 

density. Because the efficiency of grain heating by electrons decreases with increasing electron 
temperature, values higher than 2 keV have virtually no added effect on lowering the densities 
we calculate. However, lower electron temperatures can raise the inferred densities. The lowest 
observed val ue for the electro n/proton temperature ratio at a coUisionless shock is ^ .05, found 



in SN 1006 (Vinketal 



van Adelsburg et al. 



(120081) . shocks of 



120031) . In the theoretical models of [ 
greater than 1500 km s"^ have a minimum value of this ratio of 0.03. If we assume this value is 
the minimum possible value in Tycho, this leads to a minimum electron temperature of around 0.3 
keV, which is also the minir num level of electron hea ting in coUisionless shocks predicted by the 



lower hybrid wave model of 



Ghavamian et al. 



(|2007|) . Since electron temperatures this low are 
typically not found in young SNRs, this temperature can be regarded as a lower limit. The density 
inferred from a shock with kT^ = 0.3 is a factor of 1.7 higher than that inferred from a shock 
with our assumed value of 1.35 keV. Thus, while we recognize that the assumption of a constant 
value of kTg is probably not valid, it is unlikely that variations from this are large (given the fact 
that shock speeds only vary by a factor of two in the remnant), and in any case, even significant 
variations do not have a large effect on the modeling. 



3.2.3. Density Fits 

With the electron and proton temperatures estimated for each of the 19 regions, the density 
of the gas in the post-shock environment is the only remaining free parameter in our models. We 
adjust the density in each region to fit the measured 70/24 flux ratio. We show in Figure |3] the 
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70/24 flux ratio as a function of postshock density, rip, assuming constant values of kTp and kT^. 
The effect of the density on the flux ratio is quite large, allowing us to determine the density within 
a given region with relatively little uncertainty, for given values of the plasma temperature. We 
report the values of density in Table |2]and show them (normalized to the lowest density regions to 
show the magnitude of the relative density differences in the remnant) on an image of the remnant 
in Figure|2l Density uncertainties are calculated assuming a 16% error on the 70/24 flux ratio; see 
Section [U Two inferences are immediately apparent from the inferred densities. First, there are 
only three regions where the density (reported as the post-shock proton density. Hp) is higher than 
1 cm"^: two contiguous regions on the eastern limb, including the "knot g" filament seen in the 
optical, and one region in the NW. Second, aside from these three "dense" regions, there appears 
to be an overall azimuthal gradient in the densities, with the average in the NE being a factor of 
3-5 higher than in the SW. We discuss and interpret these results below, in Section HI 



4. Discussion 

The densities we determine in the various regions around the forward shock have an inverse 
relationship with the shock velocities measured in these regions, as inferred from the proper 
motion studies. This is expected, but is not an a priori constraint of the models. It is also not 
simply a result of the effect of variations in the proton temperature in these regions; as we showed 
in Section IX2l the effect of the proton temperature on the calculated IR flux ratio is fairly small, 
and certainly not enough to account for the variations seen in the density. Morphologically, 
the regions where the highest densities are measured are also the places where Ha emission 
is seen. Our models are not directly sensitive to the pre-shock ambient density, only to that in 
the post-shock environment. A standard strong shock will produce a compression ratio of four 
between the pre and post-shock gas, but this should be co nsidered a lower limit. Efficient particle 



acceleration will raise the compression ratio at the shock ([Jones & Ellison 



199lh 
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Along the eastern and northwestern portions of the remnant, the forward shock is 
encountering localized denser clumps of material. In particular, the highe st density region we find 



in the remnant corresponds to "knot g" ([Kamper & van den Bergh 



1978h . the brightest optical 



filament in the remnant. Our resu 
conclusions of 



Ghavamian et al- 



ts in these few dense knots are qualitatively consistent with the 



(|2000l) that the shock in the regions of bright Ha emission is 



encountering denser material from a neighboring HI cloud. 

However, strong Ha emission in Figure [T]is not confined just to a few dense knots located 
at the remnant's rim, but extends over the whole NE quadrant. The longest contiguous set of 
filaments in the NE is located interior to the outermost blast wave. Densities along these filaments 
are likely significantly higher than along the remnant's rim. Because of the poorly known shock 
geometry in this region of the remnant, it is not clear how the outermost blast wave related to these 
optical filaments. In one scenario, it might mark where the blast wave wraps around a low-density 
periphery of a large cloud with a high density, similar to the dense knots discussed above. 
Altematively, densities along these long optical filaments might be lower, so they would fit better 
into the framework of the overall density gradient discussed below. In this case, the magnitude 
of the overall density gradient would be underestimated with the current density measurements 
available just at the remnant's rim. Proper motion measurements of optical filaments are needed 
to resolve these ambiguities. 



4.1. Density Gradient 

The densities we determine clearly favor a low average ISM density, consistent with the lack 
of thermal X-ray emission from the blast wave. Densities in the E and NE are 3-5 times higher 
than those in the W and SW. As we will show below, hydrodynamical simulations are broadly 
consistent with this, though they suggest that the magnitude of the gradient may need to be 
somewhat higher, at around an order of magnitude. The simplest explanation for this behavior is a 
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NE-SW density gradient in the ISM. This is consistent with the hypothesis of Lee et al.1 (|2004() that 
Tycho is expanding into an I SM density gradient. It is also qualitatively con sistent with studies of 



the X-ray emitting ejecta by 



Badenes et al. 



(120061) and 



Havato et al. 



(120101) . who found that the 



ejecta are brighter in the E and NW than in the SW. Even if our models for dust emission have 
systematic errors, the fact remains that the temperature of the dust, as measured by differences in 
the IR flux ratios, varies significantly in different locations in the remnant, and variations in the 
post-shock gas density are the most plausible way to explain this. It is unlikely that our models are 
more correct on one side of the remnant than the other. Detailed multi-dimensional modeling of 
the remnant's evolution and the ionization state of the reverse- shocked ejecta is beyond the scope 
of this paper, but a full understanding of Tycho will require such work. 

A possible alternative explanation for the apparent density gradient is cosmic-ray acceleration 
that is much more efficient on one side of the remnant than the other. Our dust models are 
insensitive to this, but to explain the density differences we see purely by an increase in the 
compression ratio due to particle acceleration, the shock compression ratio in the E and NE would 
have to be > 15. While this is not beyond the re alm of possibility , it does imply that > 80% of the 



total shock energy is being put into cosmic rays (|Vink et al. 



2010l) . Also, if this were the case, then 



the ratio of the radius of the contact discontinuity (CD ) to the forward shoc k (FS) would be much 



Warren et al. 



(l2005h examined this ratio 



higher in regions of such extreme compression ratios. 
(CD/FS) around the p eriphery of Tycho, fin ding evidence for particle acceleration everywhere in 



the remnant (although 



Orlando et al 



(|2012|) have suggested that the high values of this ratio can be 
explained via hydrodynamic instabilities, without need of cosmic -ray acceleration). Furthermore, 
if such acceleration were taking place in regions of high density, we woul d expect to see a 



Warren et al. 



(l2005l) . 



correlation between the densities we measure and the value of CD/FS from 
yet we see no such correlation in our data. Even when we throw out the three highest values of 
density, those apparently coming from dense knots, we still find virtually no correlation between 
density and CD/FS (R^ = 0.14). 



-17- 



4.1.1. Gamma-Ray Emission 



Morlino & Capriolil (|2012b model the 7-ray emission from Tycho by assuming a uniform 
pre-shock density of 0.3 cm"^, finding that such a medium supports a model in which the GeV and 
TeV emission is produced by the hadronic mechanism, i.e., protons accelerated at the blast wave 
colliding with ambient pr otons to produce tt" particles, which then decay into 7-rays. However, 



Atoyan & Dermeii (|2012b find that a leptonic model of 7-rays produced via inverse-Compton 
scattering of energetic electrons off of low-energy photons fits the 7-ray data equally well. Their 
model assumes a density of 0.75 cm"^, although they note that lower densities cannot be excluded. 
Only three of the nineteen regions we model can have a pre-shock density as high as the assumed 
density in either the hadronic or the leptonic models. Clearly, further study of the 7-ray emission 
from Tycho is necessary, taking into account the density structure we find here. 



Interestingly, 



AcciarietaL 



(|2011|) reported a small offset in the TeV emission from the center 
of the remnant, in the direction of enhanced density to the NE. The TeV source is offset by 2.4' 
with a statistical uncertainty of 1.4', while the VERITAS telescope has a la PSF of 6'. However, 
even confirmation that the 7-ray morphology correlates with the density enhancements does not 
definitively select one model of 7-ray emission, as both leptonic and hadronic models could 
explain such a morphology. 7-ray emission from bremsstrahlung and 7r°-decay is enhanced by the 
higher target density, while inverse Compton scattering off low-energy IR or microwave photons 
will also be enhanced in this region. We will fully examine the implications of our density finding 
for the 7-ray emission observed in a future publication. 



4.2. An Off- Center Explosion Site 



An additional effect of an explosion into a n on-uniform ISM is that the ce nter of explosion 
will not be at the center of the resulting remnant (|Dohm-Palmer & Joneslll996|) . A recent analysis 
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by 



Kerzendorf et al 



Ruiz-Lapuente et al. 



(|2004|) within the center 



(120 12|) of six stars identified by 
of Tycho did not turn up any potential candidates for the companion star, thought to exist in the 
single-degenerate scenario of a Type la SN, leading the authors to conclude that Tycho could 
not be explained by such a scenario. However, these authors searched for potential companion 
stars only within a circle of radius 39" centered on the center of symmetry of the Chandra image. 
We report here the effects of the density gradient that we observe on the current location of the 
center of explosion with respect to the center of symmetry of the remnant. We have performed 
two-dimensio nal (2D) hydrodynamic modeli ng of Type la SN explosions with an exponential 



ejecta profile (iDwarkadas & Chevalier 



1998|) into a density gr adient, and have com pared these 



results with analytic solutions to the thin-shell approximation (|Carlton et al. 



201lh . We report 



both below. We find that the most model-independent way to infer the offset is from the observed 
velocity asymmetry, i.e., the proper motion of the shock. We assume that the observed density 
gradient is in the plane of the sky. 



4.2.1. Hydrodynamic Modeling 



in 



For our hydrodynamic modeling, we employed the same numerical methods as described 



Warren & Blondin 



(|2013|) . with the additional feature of a density gradient in the ambient 
medium. The extemal density gradient was of the form p oc e"'\ where r is a dimensionless 
distance, the same as r' in Equation (3), and l/a = rjsM is the dimensionless ISM density scale, 
discussed further in the Appendix. 

Specifically, we used the VH-1 hydrodynamics code to evolve the Euler equations for an ideal 
gas with a ratio of specific heats of 7 = 5/3 on a 2D spherical-polar grid with 300 radial zones 
by 900 angular zones, providing a spatial resolution of ^ 3.5 x 10"^ in the angular direction 
and slightly higher resolution in the radial direction. The supern ova ejecta were modeled using 



the exponential density profile of 



Dwarkadas & Chevalieii(|l998l) . We employed a moving grid to 
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track the evolving SNR as it expanded over five orders of magnitude, effectively removing any 
artifact of the initial conditions and providing sufficient time for the Rayleigh-Taylor instability of 
the contact interface between shocked ejecta and shocked ISM to reach a quasi-steady state. 

We examined various values of a to study the effects of the density gradient on asymmetries 
in the shock velocities and on the remnant shape. The values of a correspond to a range of current 
density contrasts at the current size of Tycho of factors of 5-100 from maximum to minimum. 
An important result of our simulations is that the ratio of the velocity semi-amplitude {{Y,„ax - 
Ymin)l(ymcLx + ^ min)) to the radial offset from the center of the explosion ((Rmar - Rmm)/(Rmav + 
Rmm)) is roughly constant at a value of about 2.2 ±0.1 for ages between about 300 and 700 yr. 
Our simulations predict, for different values of a and different ages, relations between the radial 
offset and the density gradient. These relations are summarized in Figure HI which shows that for 
a wide range of gradients and ages, there is a fairly tight relation between the density contrast 
across the remnant and the radial offset. One can then use either an observed density contrast or 
shock proper motions to predict the radial offset of the explosion site from the symmetry center 
for a remnant. 

For Tycho, Table [H shows that the velocity semi- amplitude of the averaged velocities from 
the radio and X-ray measurements is 0.36 (the velocities reported in that Table assume D=2.3 
kpc, but since the velocity semi-amplitude is dimensionless, it is also independent of distance). 
Dividing this by the ratio of 2.2, given above, we obtain a radial offset of 16.5% of the radius of 
the remnant, or about 40" from the geometric center of the remnant to the explosion site. The 
discrepancies between the radio and X-ray measurements are relevant here. Using the X-ray 
data alone, the velo city semi-amplitude is o nly 0.22, which leads to an offset of 25", within 



the search radius of 



Kerzendorf et al. 



(|2012|) . Using the radio data alone, we obtain a velocity 
semi-amplitude of 0.51, leading to an offset of 1'. Figure |5] shows our simulated remnant at an 
age roughly corresponding to that of Tycho, and assuming a value of a of 0.95, giving a current 
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density contrast of about an order of magnitude. This demonstrates that the remnant (in the 
relatively early evolutionary stage of Tycho) can remain remarkably round despite a significant 
external density gradient. A qualitatively similar result, for a considera bly different functional 



form of density gradient, was obtained by 



Dohm-Palmer & Jones 



( 1996 ) 



The large offset inferred from the radio proper motions alone can likely be ruled out by 
the results of our hydrodynamic and analytical modeling shown in Figure IH An offset of 1 ' 
would require a density contrast of approximately 80, well above what we measure from the 
IR data. Clearly, better proper motion measurements are necessary to resolve this discrepancy. 
We have recently been approved for new JVLA observation s, which will expand the baseline for 



(11997b to ~ 30 years. 



proper motion measurements from the 1 1 years reported in Reynoso et al. 
significantly reducing the uncertainties on the proper motions. Nonetheless, we suggest that 
searches for the progenitor compani on in Tycho be extended by a factor of two beyond the area 



searched by 



Kerzendorf et al. 



( 2012 ) 



4.2.2. Analytic Approximation 



Carlton et al. 



(|2011|) present a thin-shell approximation to a blast wave expanding into a 
uniform ISM, and here we extend their results to an ambient density gradient. There are analytic 
solutions in the thin-shell approximation for both linear and exponential density gradients (see 
the Appendix). For the linear gradient, we derive an approximate expression for the ratio of the 
dimensionless velocity semi- amplitude to the radial offset. A full derivation is presented in the 
Appendix, but the key results here are 1) the ratio from the analytic approximation is 2.4, quite 
close to the 2.2 determined from the hydrodynamic simulations above; and 2) this approximate 
expression holds even for large density contrasts inferred for Tycho where use of the exponential 
(instead of linear) density gradient is preferred. 
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4.3. Pressure 

With the density and velocity known, we can calculate the ram pressure in the shock, pov^. 
We divide the post-shock number densities by four to convert to pre-shock densities, then multiply 
that by [1.4 x (1.67 x lO"^'*)] g to get a mass density, po- We report the ram pressure for each 
region in Tabled In the east and NW knots, the pressures are a factor of 5-10 higher than 
elsewhere in the remnant, supporting the idea of a relatively recent encounter with denser clumps, 
where the pressure has not yet had time to equilibrate. Elsewhere in the remnant, pressures 
in the SW are about a factor of two lower than those in the NE, in good agreement with our 
hydrodynamic modeling. The same caveat applies here about particle acceleration increasing 
the compression ratio of the gas. If the compression ratio were higher in the NW, then the ram 
pressures of the shock could be more equal. If particle acceleration is constant in the remnant, 
only the absolute values of the pressure will be affected, but the relative differences will remain. 



4.4. Evolutionary State of Tycho 



The densities we report in Table |2] assume a distance to Tycho of 2.3 kpc and "standard" 
ISM dust grains, i.e., compact, homogeneous spheres with a mix of silicate and graphite 
grains in separate populations. If we ignore the three regions coincident with dense knots, the 
azimuthally-averaged value of the post-shock density is ~ 0.3 cm"^. Under the assumption of 
standard shock jump conditions, this leads to an average pre-shock density, Uq, of around 0.08-0.1 
cm"^. Since the IR emission is produced throughout the post-shock region where densities 
decrease from their value right behin d the blast wave, the de nsity ratio averaged over that region is 



expected to be somewhat less than 4 (W illiams et al. 



limit of 0.2-0.3 cm ^ reported by 



2011 ah. T hi s density, whil e below the upper 



Cassam-Chenai et al.. (2007 ) & 



Katsuda et al.i (^20101) . is rather 



low, given what is known about the evolutionary state of Tycho. 
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To first-order, the average density obtained for a density gradient is equal to the density being 
encountered by the forward shock in a direction perpendicular to the gradient. For Tycho, this 
would corres pond to position angles of ^ 14 and 320°. Following the scaling of hydrodynamical 



Dwarkadas & Chevalieii (|1998|) . we can write the dimensionless scaled radius, 



variables by 
r' = r/R', and the scaled time t' = t /T' zs, 



r' = 0.61 (no/0.1 cm-^f'^ {D/23 kpc) {Mej/Mch) 



-1/3 



(2) 



t' = 0.82 (no/0.1 cm-y e]',^ iM.j/MchY''^ 



(3) 



where Mg^ is the ejecta mass, Mc/, is the Chandrasekhar mass, and £51 is the explosion energy in 
units of 10^' ergs. Assuming a standard Type la explosion where Mgy = Mc/j, a pre-shock density 
of 0.08 leads t o a scaled radius of 0.57. C omparing this with the three-dimensional hydrodynamic 



simulations of 



Warren & BlondinI (|2013|) . a scaled radius of 0.57 for the forward shock implies 



a scaled time (again, using the formalism of Dwarkadas & Chevalier 1998) of only ~ 0.4. 
This is too short to explain the expansion p arameter, m (where r oc t"'), of Tycho, kn own from 



proper motion measurements to be ~ 0.53 (IRevnoso et al 



1997 



Katsuda et al 



Warren & Blondin 



20101). Recent 



(20131) show that this can be 



three-dimensional hydrodynamical simulations of [ 
achieved with a scaled time of around 1, corresponding to a scaled radius also around 1-1.1. Our 
hydro simulations confirm this; the expansion parameter and velocity semi-amplitude in Tycho 
are best matched in simulations with a scaled time of 1.05 and a scaled radius of 1.15. 

Thus, our scaled radius needs to be increased by about 70% to match hydrodynamic 
simulations of the evolutionary state of Tycho. There are three parameters in equation (2): hq, D, 
and Mej. A sub-Chandrasekhar ejecta mass would raise r', but only as M~j''^, so we view this as 
the least viable option. If distance alone were varied, a distance of 3.9 kpc would suffice. This 
would lower the density slightly (see Section 1X2] ), but since the effect on density would be small 
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and r' only goes as n^^, the overall effect of the lowered density on r' would be trivial. A distance 
of 3.9 kpc is at the high end of the range of reported distances for Tycho (see Section [T), and it 
is worth noting that such a distance combined with measured proper motions would imply that 
shock velocities exceed 6000 km s"' in some locations. 

A third, perhaps most intriguing possibility, is that the densities may be higher than we 
report in Table |2] due to alternative grain models. The assumption that dust grains are compact. 



spherical, solid bodies of homog eneous materia 



to be a physically valid model (.Shen et al. 



20081) . In 



makes for easier ca 



Williams et al. 



culatio ns, but is unlikely 



( 2011a|) . we explored the 



effect of porosity, or "fluffiness," as well as that of a non-homogeneous grain material, on the IR 
spectra produced from coUisionally-heated grains. The details of our grain model are discussed 
in that paper; we use an identical model here. For grains that are 50% vacuum, 33.5% silicate 
and 16.5% amorphous carbon, the densities we infer from the IR flux ratios are increased by 
~ 80%. This increase in the density would raise the scaled radius to 0.7. While this is not, by 
itself, enough to account for the observed evolutionary state of Tycho, it can be combined with 
a more modest distance increase to only 3.3 kpc. The resulting average density from this model 
would be 0.14-0.15, which is still below the upper limits of 0.2-0.3 discussed earlier. While this is 
suggestive that porous grains may, in fact, be ubiquitous in the ISM, current knowledge of grain 
physics is insufficient to conclusively rule on this issue. Grains of even higher porosity would 
have a larger effect on raising the d ensity, and grains exce eding 90% porosity have been suggested 



to fit the IR spectra of some stars (jLi & Greenberg 



19981). Finally, we point out that while the 



choice of grain model can affect the absolute density values one determines, the relative values 
are unchanged, and the inference of a density gradient in Tycho is unaffected. 



5. Conclusions 



Observations of Tycho in the mid-IR reveal significant color variations in the forward 
shock-heated interstellar dust. Since this dust is wanned by collisions with energetic particles in 
the postshock gas, the IR observations are a powerful diagnostic of the gas density in the remnant. 
We analyzed IR emission as a function of azimuthal angle around the periphery of the forward 
shock, finding a variation in the ratio of the 70 /im to 24 //m flux of an order of magnitude, 
implying an overall density gradient with densities higher in the NE than in the SW by a factor of 
~ 3-5. This ISM density gradient is virtually independent of distance and assumed dust grain 
composition. We also identify a few regions of significantly higher density that morphologically 
correspond to bright Ha emitting regions. These are likely localized density enhancements 
produced during the early stages of interaction with a nearby HI cloud. 

Tycho joins the growing group of Type la SNRs that are not consistent with expansion into 
a uniform ISM. We use two-dimensional hydrodynamic simulations of Type la SN explosions 
into a density gradient and find that the explosion center of the remnant must be shifted with 
respect to the geometric center of the remnant by around 10%, and possibly higher. In Tycho, 
these computational simulations agree quite well with analytic approximations. The simulations 
suggest that the magnitude of the density gradient may need to be even higher than inferred from 
IR observations, around a factor of 5-10, to explain the proper motions observed in the remnant. 
The offset of the explosion in the simulations is most sensitive to the observed proper motion 
variations of the forward shock, so improved proper motion measurements of Tycho are necessary 
for further study. Finally, simulations show that remnants expanding into a significant density 
gradient can remain remarkably round. 

The mean upstream densities in Tycho are relatively low (~ 0.1 -0.2 cm"^, depending on the 
assumed dust model), consistent with the lack of thermal X-ray emission from most places along 
the forward shock in Tycho (with the exception of the NW). The densities inferred from models 
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of standard ISM dust grains are too low to explain the evolutionary state of Tycho. Porous grains 
provide a potential resolution to this, as inferred densities from such grain models are significantly 
higher. A distance of around 3.5 kpc would also resolve this issue. Our results suggest that Tycho 
is far from a spherically symmetric, homogenous remnant, and multi-dimensional modeling is 
required for a fuller understanding. 

Facilities: Chandra, Spitzer, XMM-Newton 
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Table 1. Measurements 





OA. nTY\ PlllY 
Z-T^ LLLIL i. 1U.A. 


10 //TTl PlllY 
/ \J LLLIL i. lUA. 


70/24 Flux Ratio V, (km s"' ) T 


p i^R_c V ; 


13 


201 


1319 


6.96 ±1.1 


3660 ±366 


26.0 


31 


223 


1373 


5.97 ±0.93 


3310 ±331 


21.3 


47 


226 


1358 


5.42 ±0.84 


2360 ±236 


10.9 


63 


514 


1592 


3.50 ±0.54 


1920 ±378 


7.1 


81 


1093 


1739 


2.09 ±0.32 


2210 ±378 


9.5 


105 


279 


1394 


6.07 ±0.94 


3110±311 


18.8 


121 


315 


1441 


5.99 ±0.93 


3510 ±351 


24.0 


138 


198 


1319 


7.65 ±1.2 


3330 ±333 


21.6 


155 


181 


1361 


8.79 ±1.4 


3480 ±538 


23.6 


172 


184 


1369 


9.21 ±1.4 


3240 ±486 


20.4 


192 


103 


1281 


15.10 ±2.3 


3780 ±378 


27.8 


213 


99 


1304 


15.70 ±2.4 


4060 ±634 


32.1 


233 


94 


1339 


15.92 ±2.5 


3980 ±418 


30.9 


252 


88 


1211 


14.33 ±2.2 


3920 ±648 


29.8 


272 


96 


1248 


13.70 ±2.1 


3850 ±385 


28.8 


290 


110 


1168 


10.58 ±1.6 


3700 ±370 


26.6 


308 


346 


1499 


4.34 ±0.67 


3580 ±358 


24.9 


331 


1362 


2113 


1.98 ±0.31 


3200 ±320 


19.9 


353 


273 


1426 


5.02 ±0.78 


2380 ±238 


11.0 


Note. 


— Deg. = 


Azimuthal angle, east of north. 


All fluxes reported 


in milli- 



Janskys. Errors on fluxes are 4% for 24 /^m and 15% for 70 //m (see Section [T] 
for details). Vs = shock velocity, averaged from radio and X-ray measurements and 
assuming D=2.3 kpc. Tp = proton temperature. 



Table 2. Results 



Deg. Density (cm ^) Pressure (dyne cm ^) 



13 


0.32°:^° 


2.5l|itxl0- 


32 


0.37°:^t 


2.372|9xl0- 


47 


0.58o;4g 


1.89?j| xlQ- 


63 


1 2i-5 


2.59^:^1 xlQ- 


81 


2.1?:^ 


6.00^:2 xlQ- 


105 


o.42°:i 


2.3ll^gl xlQ- 


121 


40-50 


2.88|jf xlQ- 


138 


^■■'^0.26 


2.002;^2 xiO" 


155 


25'^-3i 


1.77^J^ xlQ- 


172 


25'^-3i 


1.53[:?^ xlQ- 


192 


12°- 15 


l.OOijgxlO- 


213 




1.06^ It xlQ- 


233 




1.021:1 X 10- 


252 


12°- 15 


LOSiif xlO- 


272 


13°i6 
u.iJo.ii 


i.i3i:^f xio- 


290 




1.52i;g xlO" 


308 


54° ''^ 


4.05^:^^ xlO" 


331 




9.58^2^^ xlO" 


353 


6^0-^1 


2.15?:^^; xlQ- 



Note. — Densities are postshock. Pressure 
calculation assumes compression ratio of 4 at the 
shock front. Both assume standard ISM dust 



grain models of lWeingartner & DrainelfcOOll) and 
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D=2.3 kpc (see text for details of dependence on 
these quantities). 
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A. Thin-Shell Solution 



Warren & Blondin 



(2013) 



For young Type la SNRs, 3-D hydrodynamical simulations by 
showed that the average radius of the contact interface between the shocked ejecta and the shocked 
ambient medium depends only weakly on details of the postshock flow. It then becomes possible 
to use a thin-shell approximation to follow expansion of the remnant. In this approximation, 
shocked ejecta and shocked ambient medium are assumed to reside in an infinitely-thin shell 
located at radius Rg at time t after the explosion and moving with velocity V^. The velocity v of the 
freely-expanding ejecta at the reverse shock is equal to Rs/t. In the exponential ejecta model of 
Dwarkadas & Chevalier (1998), the shocked ejecta mass Mej is equal to 

2 



M, 



— +2— + 2 



exp 



(Al) 



where M'"j is the total ejecta mass and Vg is the exponential velocity scale (p^y oc exp(-v/ve))- The 



shell momentum Ps is equal to the momentum of the shocked ejecta 

,2 



V 



+ ■ 



3v. 



^^^j (v/v,) +2v/v, + 2 V 
since the ambient ISM is assumed to be at rest. The shell velocity Vs is equal to the shell 

momentum divided by the total shell mass, Ps/iMej+MjsM) {Mjsm is the shocked ISM mass) 

Equation of motion for the shell is dRs/dt = V^, which can be transformed to 



(A2) 



\^dR^_ / vM,j{l+MisM/M,j) 

R, dv ~y P, 

by cha nging the independent variable from Mo v. For uniform ISM {Mism oc R^), 



(A3) 



Carlton et al. 



(|2011|) reported an exact solution to equations (IA111A31 ): 



^0 = 



1+4- + 6 - 

V V V 



1/3 



2v 



1/3 



(A4) 



where the dimensionless radius rp i s now normalized to unity for the swept-up ISM mass equal to 



M'f', as in 



(also in 



Dwarkadas & Cheva' 



Dwarkadas & Chevaliei ( 



ier J 1998 ). The dimensionless time is equal to tq (v/3'/^2ve) ^ 



1998b units). 
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Equations (|AH - IA3I) can be solved an alytically for more c omplex ISM density distributions 



than simple power laws in R discussed by 



Carlton et al. 



(|201lb . We consider here an exponential 



density distribution pisM oc exp(ar), where l/a = tjsm is the dimensionless ISM density scale, 
negative (positive) for density decreasing (increasing) with radius. For \ar\ << 1, the exponential 
density profile becomes linear, so solutions discussed below encompass both linear and 
exponential profiles. The swept-up mass is 

1 2 2 



misM = 3r 



■ + T le"'' 



ar {arf {arf J (arf 



(A5) 



in dimensionless units. By changing variables from v and R to y = {arof and x = ar m 
equations (IAlj - IA3l) . we arrive at 

x^ = Aa^misM -y = I2[{x^ -2x+2) e' -2\-y. (A6) 



A solution 



y = ^[[x^-Ax+6)e'-6-lx\ (A7) 



relates the scaled shell radius x = ar to the scaled swept-up ISM mass y = (aro)^ for ejecta 
entering the reverse shock with the same free-expansion velocity v for both exponentially- and 
uniformly-distributed ambient medium. Expansion in Taylor series at ;c = gives 

y = x' + — + j + ..., (A8) 

where the first two terms on the right-hand side provide an exact solution for the linear density 
gradient. 

Displacement Sr between solutions (IA7I) and (IA4I) . evaluated at the same time t but with 

different ejecta velocities v and vq, respectively, is given by equations 

Sr r(v)-ro(vo) Vq-v 

— = — = , VoKv) = vro(vo). (A9) 

^0 ''o(vo) V 

Generally, there is no simple explicit solution for vo and Sr, but for small displacements a Taylor 
series expansion gives 

3.(1.4..fe> 

ro 5[l + 6z(l + 3z + 4z2)]' 
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where z = Vg/vo. The dimensionless offset (Rmcix-Rmiii)/(Rmax+Rmin) between the center of 
the remnant and the true explosion center is then equal to \Sr/ro\. The dimensionless velocity 
semiamplitude (Vm^a - Vmm)/(Ktax + V^n) is also linear in a, so their ratio 

(Vn,a.-Vnun)/(V,na. + Vnnn) 3 ' ( 1 + 8z + 24f + 24z^ + 1 2^^) 

(Rmca-Rmin)/(Rmax + Rmin) [l +6z {l + 3z + 4z^)] [l +3z {l +2z + 2z^)] 

is independent of the magnitude of the density gradient. For young SNRs such as Tycho, the 
second term on the right-hand side of this equation varies only slowly with z, so at the age of 
Tycho the ratio between the dimensionless velocity semiamplitude and the dimensionless offset 
may be considered independent of age and equal to 2.4. 



Fig. 1. — Left: Spitzer 24 fim image. Center: Spitzer 70 nm image. Instrument artifacts in the 
image introduce flux uncertainties at only the 8% level (see text for details). Right: Ha image from 
2007, courtesy of P.F. Winkler. 
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Fig. 2. — The normalized density in various regions of Tycho, along with outlines of the spatial 
regions used for our analysis, superposed on a three-color IR image, with Spitzer 70 /xm emission 
in red, 24 /im emission in green, and WISE 12 /im in blue. North is up and east is to the left. The 
absolute scaling for the densities are 0.1 1 for the compact grain model and 0.2 for the porous grain 
model (see Section |4~4] and Table[2]for details). Interior numbers are azimuthal angle, in degrees. 
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Fig. 3. — The 70/24 /im flux ratio as a function of gas density, assuming electron and proton 
temperature held constant at 1.35 and 20 keV, respectively. Note that this plot is on a logarithmic 
scale on the y-axis. 
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Offset 

Fig. 4. — Offset of the explosion from the geometric center of the remnant (Rmax-Rmin)KRmwc+Rmin) 
versus density contrast in the ISM from one side of the remnant to the other. From top to bottom: 
tracks from 2D hydrodynamical simulations with values of the density gradient parameter, "a," of 
1.8, 1.6, 1.4, 1.2, 1.0, and 0.8. Diamonds show locations along each track of a dimensionless time 
of unity. The shaded gray region corresponds to dimensionless time between 0.7 and 1.3; Tycho 
should be located in this region. Plus signs mark analytic results for dimensionless time of unity, 
for values of "a" from 0.2 (lower-left) to 2.4 (upper-right). 
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Fig. 5. — Simulated remnant from a 2D hydrodynamical simulation at an age of several hundred 
years, expanding into an exponential density gradient with a = 0.95. At the time shown, the up- 
stream density is an order of magnitude higher on the left than the right. The density scale is in 
dimensionless units, and only the shocked gas is shown. The offset of the center of the explosion 
from the current center of symmetry of the blast wave is shown with an arrow, and amounts to 
roughly 10%. The outer white circle is perfectly round, and is simply shown for comparison to 
illustrate the roundness of the simulation. 
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